Orbital resonances around Black holes 
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We compute the length and timescales associated with resonant orbits in the Kerr Metric for all 
orbital and spin parameters. Resonance induced effects are potentially observable when the Event 
Horizon telescope resolves the inner structure of Sgr A*, space-based gravitational wave detectors 
record phase-shifts in the waveform during the resonant passage of a neutron star as it spirals in 
to the black hole and in the frequencies of quasi periodic oscillations for accreting black holes. The 
onset of geodesic chaos for non-Kerr spacetimes should occur at the resonance locations quantified 
here. 
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Introduction. Resonant phenomena are ubiquitous in 
multi-frequency systems and are harbingers of the onset 
of dynamical chaos [1(. In celestial mechanics, they play 
an important role in satellite dynamics. Gaps in the as- 
teroid belt and the density profile in the rings of Saturn 
[2J, |3J have in large part been sculpted by resonant in- 
teractions. The orbital motion of satellites around black 
holes is mathematically idealized as bound geodesies in 
the Kerr metric. Unlike in Newtonian gravity where or- 
bits are characterized by a single rotational frequency u;^ , 
geodesic motion in the Kerr spacetime is described by 
three frequencies [J]. The two libration type frequencies 
oj r and ujg corresponding to motion in the radial and Ion- 




Event Horizon 

FIG. 1: (Color online) The location of low order resonances around 
a black hole superimposed on an embedding diagram. The line 
width of each resonance is inversely proportional to the order of 
the resonance to give an indication of the relative importance of a 
particular resonance. 



gitudinal directions augment the azimuthal frequency ui$. 
The existence of these libration type frequencies gives rise 
to the resonant phenomena considered here. 

When exploring dynamics in an astrophysical environ- 
ment such as near Sgr A* at the Galactic center a number 
of corrections to the vacuum Kerr Hamiltonian Hk must 
be taken into account. The presence of an accretion disk 
|5|, other sources of matter, structural deviations of the 
central black hole away from the Kerr metric J6|-4s|l , the 
influence of modified gravity, or the mass and spin of the 
satellite [9, [l0( will all affect its orbital motion. Regard- 
less of the nature of the perturbation, the Kolmogorov- 
Arnold-Moser (KAM) theorem states that the perturbed 
dynamics will be a smooth distortion of Kerr geodesies 
provided the frequencies of the motion in "Hk are suffi- 
ciently irrational as quantified by the criterion ll|, [12 1 
\mu) r — nujg\ > K (e) / (n + m) 3 . The factor K(e) here ap- 
proaches zero as the perturbation vanishes. The notable 
exception to this theorem is low order (n + m value) reso- 
nant orbits whose frequencies occur in the rational ratios 
bj r /ijjg = n/m = 1/2, 1/3, 2/3, • • • . For these orbits the 
possibility of dramatic deviations from Kerr dynamics ex- 
ists. Since the predictions of the KAM theorem depend 
on %k only we expect a measurable imprint of Kerr's res- 
onant structure in any astrophysical environment. The 
relative locations of low order resonances is shown in the 
embedding diagram in Fig. Q] and the characteristic time 
and length scales associated with each resonance tabu- 
lated in Table |U For Sgr A*, the typical timescales for 
the low-order resonances are ~ lhr and the resonances 
occur at distances of ~ 50^as from the black hole. 

Within the next ten years radio telescopes will attain 
sufficient sensitivity to resolve length scales typical of 
resonant phenomena at the center of our galaxy. The 
event horizon telescope is one such observational effort 
currently under development [13| using the VLBI net- 
work. As a neutron star or stellar mass black hole spirals 
into a super-massive black hole, it will sample all the res- 



onant bands depicted in Fig. [I] Future space based grav- 
itational wave detectors may observe resonance-induced 
shifts in the phasing of the emitted gravitational waves 

l4 Il5j . X-ray, optical and infrared telescopes do not 
have the resolving power to image Sgr A* directly but can 
potentially record flux variations from this region that 
may display timescales characteristic of resonant events 
[19 ] . Quasi-periodic oscillations (QPOs) observed in the 
X-ray spectra of several black hole candidates that ex- 
hibit peaks at frequencies in a low integer ratio could 
also potentially be associated with the orbital resonances 

171 [l8| . To aid the identification of astrophysical phe- 
nomena that originate from orbital resonances we fully 
characterize the region of parameter space where reso- 
nant effects occur. We present a number of easily eval- 
uated formulae demonstrating the spin and eccentricity 
dependence of resonances and build an intuitive under- 
standing for the inclination dependence. 

The resonance condition. Geodesic motion in the 
Kerr spacetime with spin parameter a, is integrable. 
The energy E, azimuthal angular momentum L z and 
Carter constant Q fully specify the trajectory of a par- 
ticle with rest mass /i [l9j |. An equivalent description 
of the trajectory can be found using Kcplcr-typc vari- 
ables that are directly related to the orbit's geometry [4J: 
its semi-latus rectum, eccentricity, and the cosine of the 
maximum orbital inclination, {p, e,z_}. For a generic 
bound orbit expressed in Boyer-Lindquist coordinates 
(t, r, 8, <fi), the radial motion oscillates between the apas- 
tron, ri = p/(l — e), and the periastron, r 2 = p/(l + e), 
with a frequency u r . The longitudinal motion oscillates 
about the equatorial plane with a frequency uig, sampling 
the angles —6* < 8 < 8*. we define z_ = cos 6**. 

Resonances occur for parameter values on a two- 
dimensional surface in {p, e, Z- } space determined by 
the resonance condition 
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where the functions R and can be factored as 

R = -/3 2 (r - n)(r - r 2 ){r - r 3 )(r - r 4 ) 
G = a 2 (3 2 (z 2 -z 2 „)(z 2 -z 2 + ). 



(1) 



(2) 
(3) 



Here, (3 2 = (/i 2 — E 2 ) and the roots obey r\ > r 2 > 
r,3 > r4 and Z-\- > Z-. Evaluating the right-hand side 
of (QJ and extracting the physics of the resonant surfaces 
is complicated by the fact that the roots r 3 , r^ and z+ are 
implicit functions of {p, e, z_}. The resonance condition 
([1]) can be expressed in its most symmetric form using 
Carlson's integrals [21| which allows us to obtain several 
useful analytic results in various limits [20|. A rapidly 
convergent semi-analytical scheme for finding these sur- 
faces in general can be constructed using these integrals 

m. 



Eccentricity dependence 




FIG. 2: (Color online) The location of the 2/3 resonance in 
{p, e,z_} parameter space. The arch shape is typical for all reso- 
nances at fixed spin and depends very weakly on eccentricity. A 
maximum value of z^_ = 1 is reached at p = p po i a r ~ P* = 10.8. 
For a given e, the maximum (retrograde) and minimum (prograde) 
values of p occur on the equatorial plane z_ = and are indicated 
in blue (right) and red (left) respectively. 



Features of resonance surfaces. The 2/3 resonance sur- 
face in {p,e,Z-} space for a maximally spinning black 
hole is illustrated in Fig. [2] For a given spin, all res- 
onant surfaces display the same qualitative eccentricity 
and inclination dependence. The resonance surface has 
the shape of an inverted 'U' arch that is weakly depen- 
dent on eccentricity and attains a maximum inclination 
of zt = 1 at p = pp i ar , indicated by a black line in Fig. [2] 
For smaller inclination, z 2 _ < 1 and fixed eccentricity, the 
two possible values of p on the resonant surface corre- 
spond to prograde, p+ < p po iar, closer to the black hole 
and retrograde, p- > p po iar, resonant orbits respectively. 
The p± subscript identifies sgn(aL z ) = ±1. As z- de- 
creases the distance (p_ — p+) monotonically increases to 
its maximum value on the equatorial plane. 

Since the basic features of a resonance for a given spin 
parameter a depend weakly on eccentricity, we will now 
choose a representative eccentricity and study the spin 
dependence. Fig. [3] shows the 2/3 resonance surface for 
an eccentricity of e = 0.5 as a function of black hole spin 
and p. Here we see that the arch-width exhibits a strong 
spin dependence, with its inverted 'U' profile pinching off 
to a single column 'I' profile at p = p po iar when a — > 0. 
This indicates that resonances in Schwarzschild space- 
time are independent of inclination because the longitu- 
dinal frequency is just the same as the (^-frequency in 
this case. As the spin increases the opening angle of 
the arch increases to a maximum arch width for a = 1. 
The result is a '^'-shaped footprint of the arch in the 
p. a plane. As inclination increases, the prograde and 



Spin dependence 
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32.7 


5.10 


1/2 
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40.9 


50.3 


3.31 


1/3 
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110.2 


34.5 
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4.27 
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10.8 


223.0 


55.2 


78.9 


2.11 


1/4 


6.4 


101.7 


32.7 


36.0 


4.63 


3/4 


13.7 


319.1 


70.1 


112.9 


1.48 


1/5 


6.3 


98.2 


31.9 


34.7 


4.80 


2/5 


7.1 


119.9 


36.5 


42.4 


3.93 


3/5 


9.4 


180.4 


47.9 


63.8 


2.61 


4/5 


16.7 


427.5 


85.1 


151.3 


1.10 



1 4 



TABLE I: Time and length scales associated with low-order res- 
onances. The values are for the e = a = z— = vertexes seen in 
Fig. [4] both in dimcnsionlcss units scaled by the black hole's mass 
M and in physical units for the special case of the Galactic center, 
Sgr A*. The values are calculated from p* = 6/[l — (n/m) 2 ] and 
T = 2vrp* 3 / 2 . 



FIG. 3: (Color online) Spin dependence for the 2/3 resonance 
with fixed e = 0.5. Observe that as a — > the arch pinches off 
to a line at p = p po i ar , indicating that resonances in this limit are 
independent of inclination. The maximum arch width occurs for a 
maximally spinning black hole a = 1. 
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FIG. 4: (Color online) The location of low order resonances (m < 
7) for e = z— = as a function of a and p. For a = the prograde 
and retrograde branches are degenerate at p = p* . Each vertex is 
labeled by the rational ratio n/m and darker colors indicate lower 
order resonances. As the spin increases the retrograde branch leans 
right (copper tinge) and the prograde branch leans left (blue tinge). 



retrograde branches of the arch approach p po iar and the 
'V' profile narrows from its largest opening angle on the 
equatorial plane to a line for z~ = 1. 

A key result derived in [20J is that an exact analytic 
solution for the 'V' profile of circular equatorial orbits can 
be found. Since the 'U' profile is only weakly dependent 
on eccentricity this allows us to benchmark the location 
of the n/m resonance for a black hole of arbitrary spin. 
When e = z_ = 0, Eq. (P) is equivalent to 

[p(p-p*) - a 2 (p* - 3)] 2 - 4a 2 p(P* - 2) 2 = 0. (4) 
where p* identifies the resonance we are considering and 



is related to the ratio n/m by 



6 



1 - [n/mf 



(5) 



For non spinning black holes, p = p* is solution to Eq. (U) 
and determines the position of the '/' column in the 'U'- 
' V transition for circular orbits. The value of p* sets 
the general mean radius in physical space about which 
all the interesting features associated with the n/m res- 
onance occur. Numerical values of p* for several low or- 
der resonances are given in Table. |U For spinning black 
holes, the largest two roots of Eq.(j4|) yield the 'V pro- 
file on the equatorial plane. This profile is plotted for 
several low-order resonances in Fig. [4] The maximum 
splitting of the retrograde and prograde branches occurs 
when a = 1 in which case the relevant roots of Eq. ((4|) are 
p± = p* — 1 T 2-y/p* — 2. For low spin values, the series 
expansion 



p T =p* ± 



2a{p* - 2) 



a 2 (p* 2 



5p* + 8) 
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is useful for making astrophysical estimates. The weak 
eccentricity dependence of the 'V profile for a = is 



£-1 

p* 



e 4 (4p* - 17) 



A(p* - 6) 64(p* 



+ 0(e 6 ). (7) 



Observe that as the resonant surfaces approach the in- 
nermost stable circular orbit (p — 6) the effects of eccen- 
tricity become increasingly important. 

Astrophysical implications. Generally resonances can 
have either a capturing or a destabilizing effect on parti- 
cles that enter their region of influence. We expect Kerr 
resonances and the 2/3 resonance in particular to have 
a stabilizing influence under perturbations [201] . For par- 
ticles that are light enough, entering a resonance zone 
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FIG. 5: (Color online) The 2/3 resonant surface for a = 1 pro- 
jected onto E, L z coordinates. Lines indicate constant values of 
z£. = {0,sin j, l/\/2, cos 2., 1}, and e = {0, |, ^, |, 1} correspond- 
ing to the colored lines in Fig. [2] The L z = line corresponds to 
the 2_ = 1 line in Fig. [2] the sharp cusp for L z > to the prograde 
e = 2_ = point on the arch (left front), the lower corner point 
for L z < to the arch's retrograde e = z_ = point (right front) 
and the E = 1 line to e = 1 for the arch. The large asymmetry 
across the L z = plane is as result of the high spin value. 



can temporarily halt the radiation-reaction driven orbital 
evolution because of an enhanced exchange of energy and 
angular momentum from the black hole to the particle 
orbit. If a particle becomes captured by a resonance its 
orbital parameters are expected to change within the res- 
onant surface [22J and for generic perturbations the orbit 
should simply evolve to the lowest energy state [23|. In 
Fig. [5]we show the orbital energy and azimuthal angular 
momentum associated with the resonance surface that 
was depicted in {p, e, Z-} space in Fig. [2] The lowest 
orbital energy state for the given resonance occurs in the 
lower right-hand corner of Fig [5] corresponding to pro- 
grade circular equatorial orbits. The migration of reso- 
nantly capture particles to circular equatorial configura- 
tions could result in a cohesive resonant structure where 
the ejection of particles due to collisional interactions is 
minimal. This should result in a distinctive imprint of 
density inhomogeneities on any thin disk surrounding a 
black hole similar to the structure imprinted on Saturn's 
rings [2J, |25j . Unlike the rings of Saturn whose structure 
is largely stationary, the black hole rings would be dy- 
namical because the effect of radiation reaction dissipa- 
tion will alter the resonant structures and enable escape 
from resonance [22J, [26j . Once a trapped over-density be- 
comes sufficiently large the radiation reaction force will 
dominate over the resonance's trapping potential and the 
ring will break, depositing an over-density of matter on 
the next ring in the disk. Thus in one related catastrophic 
event the whole ring structure will re-adjust. The radi- 
ation emitted in the process is likely to be modulated 
with the characteristic frequencies associated with the 
collapsing resonance bands. The X-ray spectrum of a 
black hole candidate shows QPOs at pairs of frequencies 



in a 3:5 ratio in addition to the 2:3 ratio also observed 
in other systems |l7| . Observing a 3 : 5 frequency ratio 
is unexpected; from a dynamical systems perspective it 
is expected to be weaker than the 3:4 resonance. A pos- 
sible explanation for its occurrence is that the physical 
location for the 3/5 resonance (~ 9.4M) occurs just in- 
side the 2/3 resonance (~ 10.8 M) for all spin values as 
shown in Table U and Fig. |U Matter from a disruption 
event of the 2/3 resonance structure could collide with 
even a tenuous over-density of matter at the 3/5 reso- 
nance location due to their close physical proximity and 
stimulate photon emission. 

The results presented in this paper may provide a ro- 
bust method of determining the black hole's spin given 
observational evidence from more than one resonance. 
Recent monitoring of Sgr A* with the 1.3mm VLBI 
showed time- variable structures on scales ~ 4 R s [16|, [27 1 . 
The physical origin of this structure is not yet clear, but 
the scale is similar to that of the low-order resonances 
given in Table |U For argument's sake suppose that the 
origin of the structure at ~ 4i?s = 8M is due to the 2/3 
resonance that is displaced from its non-spinning position 
listed in Table U due to Sgr A*'s spin. On astrophysical 
grounds the 2/3 resonance is likely to have the great- 
est probability of being directly observable [20|. Using 
Eq. © the spin displacement of the prograde 2/3 res- 
onant surface is p+/M = 10.8 — 5.36a, and comparing 
to the observed structure indicates that Sgr A* has spin 
a = 0.5. The plausibility of identifying this structure 
with the 2/3 resonance could be confirmed if characteris- 
tic timescales of slightly less than an hour are associated 
with the variability and a 2:3 ratio in observed frequen- 
cies is discovered. Note that once the spin is determined 
Eq. ((6J predicts the location of the other resonances. As 
the sensitivity of the VLBI measurements increase and 
more of the horizon scales are resolved the observation of 
further resonances provides an independent check on the 
spin determined using the 2/3 resonance. 

Long-term monitoring of time of arrival signals from 
a pulsar with orbital period of a few months with the 
Square Kilometer array could determine the mass, spin 
and quadrupole moment of Sgr A* to a precision of 
< 10 -2 [28| |. providing an exciting prospect for a defini- 
tive test of the no hair theorems. A corollary of the 
results in this paper is that orbits with periods of or- 
der months are sufficiently far from the low-order reso- 
nances that the KAM theorem guarantees the region to 
be effectively free of stochastic motion. Frequency drifts 
computed using perturbative methods based on averag- 
ing as done in [29( thus accurately describe the physical 
system. Tracking the trajectory of a pulsar in the region 
50R S < r < 1000-R S should build up an accurate map of 
the central object's gravitational potential. 

From Table H] we observe that future space-based 
gravitational wave detectors sensitive to frequencies 
~ 10~ 4 — 10 _1 Hz will directly probe resonant dynamics. 



This is an exciting possibility but it underscores the ne- 
cessity of carefully modeling and incorporating resonant 
effects in the search templates. The resonance induced 
phase shift could potentially make careful parameter es- 
timation difficult. If the central object is a non-Kerr 
Black hole the possible onset of geodesic chaos will occur 
first in these regions and further complicate the analysis. 
Further numerical investigation to fully quantify these 
effects for all E and L z values, given in Fig. [5] and [2flJ, 
associated with resonance surfaces is required. 

Conclusion. We have explored the basic properties of 
resonant surfaces associated with radial and longitudinal 
motion around a Kerr black hole. A few simple expres- 
sions to quantify resonant effects in astrophysical sys- 
tems with strong gravitational fields have been provided. 
We have suggested a resonance-based method for deter- 
mining the spins of black holes. Observations of QPOs, 
gravitational wave emission from resonant transits and 
radio maps of Sgr A* at event horizon scales will in the 
near future provide a powerful observational toolkit for 
probing resonant phenomena. 
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